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Abstract 

Using renormalization-group methods, differential equations can be obtained for the all- 
orders summation of leading and subsequent non-leading logarithmic corrections to QCD 
perturbative series for a number of processes and correlation functions. For a QCD per- 
turbative series known to four orders, such as the e + e~ annihilation cross-section, explicit 
solutions to these equations are obtained for the summation to all orders in a s of the leading 
set and the subsequent two non-leading sets of logarithms. Such summations are shown for 
a number of processes to lead to a substantial reduction in sensitivity to the renormalization 
scale parameter. Surprisingly, such summations are also shown to lower the infrared singu- 
larity within the perturbative expression for the e + e~ annihilation cross-section to coincide 
with the Landau pole of the naive one-loop running QCD couplant. 
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1 Optimal RG Improvement of F(B — > X u i vj) 

Optimal renormalization-group (RG) improvement of a perturbative series to a given 
order in the expansion couplant is the idea of including within that series all higher-order 
contributions that can be extracted by renormalization-group methods [1]. We call such 
terms, which involve leading and successive logarithms of the renormalization scale /x, 
RG-accessible. Techniques have been developed to obtain closed-form summations of 
such RG-accessible contributions to all orders in the perturbative expansion parameter, 
and such summations have been shown to lead to "optimally RG-improved" expressions 
for perturbative quantities that have significantly diminished dependence on fj, [2, 3]. 

For example, leading and subleading perturbative QCD corrections to the inclusive 
semi-leptonic B — > X u t~vi decay rate, which in tree-order is purely a charged-current 
weak interaction process, are given by a QCD series [4] 
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If one substitutes eq. (1) into eq. (2), one obtains a decay rate that decreases 
monotonically with increasing //, raising the question as to which value of fi is most 
appropriate for comparing the calculation (2) to experiment. Clearly, such dependence 
on the unphysical parameter fx is an embarassment; indeed the renormalization group 
equation for the series S 
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is nothing more than a chain rule expression for the requirement that the physically 
measurable decay rate be impervious to changes in the renormalization scale parameter 
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The residual ii-dependence of the decay rate obtained from the series (1) is necessarily a 
consequence of the truncation of that series, as well as the relatively large value of the 
expansion constant a s (fi)/n. 

In fact, the series (1) may be expressed as a double summation over powers of loga- 
rithms and the expansion parameter, i.e., in the following form: 
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where 

x = a s (fj) /tt, L = \og(fi 2 /m 2 b (fi)). (6) 

The first few constants of this series, i.e., the set {T 0i0 (= 1), T lj0 , T 11? T 2i0 , T 2j i, T 2)2 }, 
are given by eq. (1). However, all higher-order constants of the form T n ra , T n>n _i and 
T n>n -2 can be obtained via eq. (3), and, hence, are RG-accessible. In terms of the new 
variables x and L, the RG-equation (3) may be expressed as 
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If we substitute the series (5) into the RG-equation (7), as well as the known QCD series 
expansions of the RG-functions 
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(3(x) = -Y.PnX n+2 , lm {x)=-Y.lnX n+ \ (8) 
n=0 n=0 

we find for any integer p that the aggregate coefficients of x p L p ~ 1 , x p L p ~ 2 and x p L p ~ 3 
on the right hand side of eq. (7) necessarily vanish: 

x p L p ~ l : 0=pT P:P -(3 T p _ 1 , p _ 1 (p-l)-5 l0 T p - 1 , p - 1 (9) 

x p L p ~ 2 : = (p-l)T P:P _ 1 + 2 l0 (p-l)T p _ 1:P _ 1 -f3 (p-l)T p _ 1:P _ 2 

- Pi(p ~ 2)Tp_ 2jP „ 2 - 57oT p _ liP _ 2 - 57iT p _ 2jP _ 2 (10) 

x p L p-z. o = (p-2)T^_ 2 + 27o(p-2)T p _ l!P _ 2 + 2 7l (p-2)T p _ 2iP _ 2 

- p (p ~ l)T p _i iP _3 - P ± (p - 2)T p _ 2jP _3 - p 2 (p - 3)T p _ 3iP _3 

- 57 T p _i iP _3 - 57iTp_ 2i p_ 3 - 57 2 Tp_ 3i p_ 3 . (11) 

Given knowledge of T 00 (= 1), one can calculate any coefficient T pp through successive 
applications of eq. (9). Indeed the eq. (1) values Ti i = 5 and T 22 = 415/24 follow 
from just two successive iterations of (9) using the rif — 5 QCD values 70 = 1, (3q — 
23/12. Similarly, knowledge of all T P:P coefficients, as obtained via (9), plus knowledge 
of T 10 = 4.25360 [eq. (1)] is sufficient via successive applications of (10) to determine 
all coefficients T PjP _i. Finally, knowledge of all coefficients T P:P T p , p _i plus the single 
coefficient T 2 = 26.7848 is sufficient via successive applications of (11) to determine all 
coefficients T PiP _ 2 , since the set of MS RG-function coefficients Pa,Pi,Pi,lo{= l),7i. 
and 7 2 } have all been calculated [5]. 

Since we now see that all coefficients T PtP) T PtP _i and T PiP _ 2 are RG-accessible, it 
makes sense to restructure the double- summation series (4) in the form 

00 00 00 00 

S[x, L] = Y, T p , p x p L p + ]T T p ^ lX p L p - 1 + ]T T PtP _ 2 x p L p - 2 + ^ T PtP _ 3 x p L p - 3 + ... , 

p=0 p=l p=2 p=3 

(12) 
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since the first three terms above are completely determined by eqs. (9), (10) and (11). 
We express (12) in the more compact form 



S[x,L] = J2x n 



n=0 lp=n 
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= Y^^SnW (13) 
n=0 

and note that S (xL), Si(xL) and S 2 (xL) all correspond to RG-accessible functions, 
based upon the information given in (1). Indeed, the program of optimal RG-improvement 
is nothing more than the explicit closed-form evaluation of these functions, and their 
subsequent incorporation into the calculated decay rate. 

To evaluate the summation S (u), as defined by (13), we simply multiply eq. (9) by 
u p ~ l and sum from p — 1 to infinity: 

oo oo 

p=i p=i 

oo 

= (1-^)^- 570^0. (14) 

We note from the definition (13) of the series S n (xL) that 

S n (0) = T n , . (15) 

The solution of the differential equation (14) with initial condition So(0) = T 0j o = 1 is 

S (u) = (l-p ur 5 ^ fS °. (16) 

A similar procedure is employed to find Si, and S 2 . If we multiply eq. (10) by u p ~ 2 
and then sum from p = 2 to oo, we find after a little algebra that 

(1 - M^- - (ft + 5 l0 )Si = 5 7 iS - (2 7o - M^. (17) 

Substituting the solution (16) into the right hand side of (17) and noting that Si(0) = Ti ;0 , 
we find that 

5(7oA/A) - 7i)/A) 



Si(u) = 



(1 - p u) 5 ™/P° 



Ti,o - 5(70/3x7/30 - 7 i)/A) + [57o(27o - Pi/Po)/(h] Ml ~ M hR) 

(1 _ /5oW )(/3o+57o)//3o • ^ 

Similarly, we can multiply eq. (11) by u p ~ 3 and then sum from p = 3 to oo to obtain the 
differential equation 

_ (2A) + 5 7o )S 2 

= ((3 lU - 2 7o )^ + (Pi + 571)51 + (P 2 u - 2 7 i)^ + 5725o(19) 
whose solution is given by eq. (2.28) of ref. [2]. 
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2 Order-by-Order Elimination of Renormalization 
Scale Dependence 

If one substitutes solutions for So(xL), S\(xL) and S^xL) into eq. (13), one obtains 
the following optimally RG-improved version of the series S [2]: 



S (xL) + xS ± (xL) + x 2 S 2 (xL) 
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with 



1 — {3qxL 



(20) 



(21) 



and with x and L given by eq. (6). 

When multiplied by ml(fi), this expression has the remarkable property of being almost 
entirely independent of /i. Figure 1 of ref. [2] displays a head to head comparison of the 
fj, dependence of eq. (20), and the same expression with S given by the known terms of 
eq. (1). For the latter case, [mj(/i)] 5 5 is seen to decrease from xs 2500 GeV 5 to 1500 
GeV 5 as // decreases from 1.5 GeV to 9.0 GeV. For eq. (20), however, the quantity 
[m p (/j,)] 5 S is seen to be 1816 + 6 GeV 5 over the same range of fx, effectively removing all 
yu-dependence from the optimally RG-improved two-loop calculation. 

Such elimination of renormalization-scale dependence via optimal RG-improvement is 
also upheld for a number of perturbative expressions, including QCD corrections to the 
inclusive semileptonic decay of B-mesons to charmed states (B — > X c t~vi), QCD correc- 
tions to Higgs boson decays, the perturbative portion of the QCD static potential function, 
the (Standard-Model) Higgs-mediated WW — > ZZ cross section at very high energies, 
and QCD sum-rule scalar- and vector-current correlation functions [2]. This last example 
is of particular relevance for QCD corrections to the benchmark electromagnetic cross- 
section ratio R(s) = a(e + e" — > hadrons)/(j(e + e~ — > /i + /i~). Such QCD corrections 
are proportional to the imaginary part of the vector-current correlation function series, a 
series which is fully known to three subleading orders in a s [6]. For five active flavours, 
we have 

S = 3R(s)/U 

= l + x+ (l.40924 + x 2 + (-12.8046 + 7.81875L + j^ 2 ) x 3 + ... 

(22) 

where x = a s (fx)/n, as before, and where L is now the logarithm 

L = log(/i 2 /s). (23) 
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Note that dependence on the physical scale s resides entirely in the logarithm, and 
that the all-orders series (22) for S, a measurable quantity, is necessarily impervious to 
changes in /i. However, progressive truncations of (22) introduce progressively larger 
amounts of renormalization scale dependence. For example, if the series S is truncated 
after all its known terms, as listed in eq. (22), we find for y^s = 15 GeV that to order x 3 , 
S increases modestly from 1.0525 to 1.0540 as /j, increases from 7.5 to 30 GeV. 1 Had 
we truncated the series (22) following its 0(x 2 ) term, we find that such a truncation of 
S now decreases from 1.056 to 1.053 over the same range of \i, doubling the magnitude 
of /i-dependence evident over this range. Finally, if we consider only the lowest order 
correction to unity (S — 1 +x(fi)), we find that S decreases from 1.061 to 1.045 as fi 
increases from 7.5 GeV to 30 GeV. 

Optimal RG-improvement of the known terms of the series (22) has been shown by 
the same methods delineated above to lead to the following expression [2]: 



S = 1 + x/w + x 2 [1.49024 - 1.26087 log w]/w 2 



+ x 3 



0.115003^- 12.9196 -5.14353 \ogw + 1.58979 log 2 w /w 3 + ... (24) 



where w is given by the definition (21), but with L now given by eq. (23). Eq. (24) is, of 
course, really a restructured version of the same infinite series as eq. (21), and similarly 
must be independent of fj, when taken to all orders. However, for the series (24) such 
imperviousness to changes in renormalization scale is now evident on an order-by-order 
basis. Truncation of the series (24) after its first nonleading term (i.e., S = 1 + x/w) 
still provides an expression that exhibits less variation with /x than all four known terms 
of the series (22). As /i increases from 7.5 to 30 GeV, we find for y/s = 15 GeV that 
1 +x/w decreases from 1.0524 to 1.0516. Similarly, truncation of the series (24) after 
its third term leads to a slow decrease from 1.0557 to 1.0553, and retention of all four 
known terms leads to an almost flat value (1.05372 ± 0.00004) over the same 7.5 GeV - 
30 GeV spread in ji. 

Consequently, we see that the program of optimal RG-improvement, as described 
above, is seen to yield order-by- order perturbation-theory predictions which are almost 
entirely decoupled from the particular choice of renormalization scale. 



3 Lowering the Infrared Bound on Perturbative 
Approximations to R(s) 

The optimally RG-improved series (24) is term-by-term singular when w — [1 — (3qx(h) 
log(/i 2 /s)] is zero. Since s is the external momentum scale characterising the physical 
e + e~ annihilation process, we see that the use of (24) is possible only if [3] 

a ^f exp (-^y ■ (25> 

1 In all estimates presented here, a s (fi) is assumed to evolve via its known 4-loop order /3-function 
from a s {M z ) = 0.118 [7]. 
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It is particularly curious that this bound on s corresponds to the infrared bound on the 
naive one-loop (1L) running couplant (x 1L = (a s (fi))i L /ir) 



^ = ~(3oxl L (26) 

whose solution 

a ' M = »VAfj (27) 

can be inverted as follows to express the 1L infrared cut-off in terms of some reference 
value of a s (n): 

A ^M™)- (28) 

Consequently, for a given choice of \i for which a s (/i) is known (e.g. the value a s (m T ) 
extracted from r-decay experiments), we see that each of the terms in the optimally 
RG-improved series (24) diverges as s approaches A 2 L from above. 

The idea that QCD perturbative series break down in the infrared is hardly new, but the 
location of this breakdown is usually identified with an IR-divergence in the all- available- 
orders evolution of a s (fi), not the naive 1L Landau pole of eq. (28). To consider the 
IR boundary of QCD corrections to e + e~ annihilation, we find for three active flavours 
that QCD corrections to the e + e~ annihilation cross-section are, as before, obtained from 
the perturbative series within the imaginary part of the QCD vector current correlation 
function [6]: 



R(s)/2 = l + x+ ^1.63982 + ^x 2 



+ (-10.2839 + 11. 3792L + ^L 2 ) x 3 

+ ... (29) 

where x = a s (fi) and where L = log(yU 2 /s). The standard phenomenological approach 
to this series is to first recognize its all-orders invariance under changes in ji, and then 
to assume such invariance applies to truncation of the series after its four known terms. 
This (seldom stated) assumption [8] motivates the choice /i 2 = s (i.e. L = 0) leading to 
the usual rif = 3 expression [7], 

R(s) = 2 [l + x(y/s) + 1.63982 x 2 (v / 2) - 10.2839 x 3 (y/s) + ...] . (30) 

Such an expansion necessarily falls apart in the infrared when x(y/s) becomes large. 
Indeed, the large coefficient of x 3 (y / i) in (30) manifests itself in a sharp drop in R(s) at 
v^ = 650 MeV [3]. 

It is interesting to compare the known terms in eq. (30) to the optimally RG-improved 
version of the known terms in eq. (29) [2, 3]: 

R(s) = s[l + x(fj,)/w(fj l ,s) 
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where 



+ x 2 (li) (l.63982 - — log (w(li, s))J /w 2 (/i, s) 

+ x 3 (^) (-1.31057^(^,5) -8.97333 

- 8.99096 log (w(/i, s)) + 3.16049 log 2 (w(/i, s))) /w 3 (/i, s) 
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w(/i, s) = 1 - -x(/i) log(^ 2 /s). 



(31) 



(32) 



We first note that x(li) and x(a/s) occurring in eqs. (29), (30), (31) and (32) are 
evolved through all known orders of the /9-function (8). Consequently, we are free to 
assign to eq. (31) an rif = 3 empirical value for fx (fx = m T or, alternatively, fx = 
1 GeV) safely outside the infrared region. To facilitate comparison of eqs. (30) and 
(31), we will assume that x(y/s) devolves via the full /3-function from this same initial 
choice of fx until, for a sufficiently small value of s, x(y/s) becomes infinite. The point 
here is that for nj = 3, the first four known terms of the /5-function are all same-sign: 
(Pox 2 + P±x 3 + p 2 x 4 + fax 5 ) > (3 x 2 . Thus for a given value of x, the full /3-function 
is more negative than the one-loop /^-function of eq. (26). Since the evolution of both 
equations is referenced to the same initial value fx, the all-orders couplant x(y/s) will 
diverge at a value of s that is larger than K\ L , the Landau pole of eq. (26). Hence the 
series (31) will probe more deeply into the infrared than the series (30) for R(s). 

As an example, consider the running couplant x(y/s) obtained via a two-loop (3- 
function f3(x) = —(3qX 2 — (3ix 3 . Solution of the differential equation s dx/ds = f3(x) 
with the initial value x([x) yields the exact constraint 




1 



1 



x(fx) 



ft, 

( r\ v IT lo S 

X{y/S) (Jq 



X(fl)[x(y/s) + Po/Pl 
X(y/S)[x(n) +P /Pl 



(33) 



The two-loop Landau pole s 2 l occurs when x( v /s 2 T) — > oo, i.e., when 

/9i//3g 



s 2L = At exp 



1 



1 + 



Pix(fx) 



A 2 



1 + 



Pix(fx) 



Pi I Pi 



(34) 



Since f3 , Pi and x(fx) are all positive, s 2 l > A 2 L . Note that p 2 and p 3 [9] persist in 
being positive. Consequently, for a given initial value x(fx), the singularity in eq. (31) 
that occurs at w(fx, s) = (i.e., at s = A 2 L ), continues to precede the Landau singularity 
of x(y/s) characterizing eq. (30). Thus, the optimally RG-improved eq. (31) extends the 
applicability of perturbative QCD to lower values of s than in the conventional eq. (30) 
approach to R(s), as is explicitly shown in Fig. 2 of ref. [3]. 

Finally, we note that one must distinguish between the infrared limitations on the 
domain of perturbative approximations to R(s), and any such limitations on R(s) itself. 
For example, each term within the toy series S^Lo [ — x /( s ~ A 2 )]" diverges at s = A 2 , 
but the function (s — A 2 ) /(s — A 2 + x) from which this series is extracted, is clearly finite 
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at s = A 2 . Similarly, the all-orders function R(s), as opposed to truncations of its series 
representations, may indeed proceed smoothly to a finite limit as s — > [10]. If such is 
the case, the best one can hope for in a perturbative series representation of R(s) is the 
deepest possible penetration of that series into the low-s region. 
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